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Abstract 



O 
O 

o 

(N 

D ' The exclusive rare radiative B meson decay to the orbitahy excited tensor 



X|(1430) meson is investigated in the framework of the relativistic quark 



l^ . model based on the quasipotential approach in quantum field theory. The 

■ calculated branching ratio BR{B -^ /s:|(1430)7) = (1.7 ± 0.6) x 10"^ as weh 

as the ratio BR{B -^ K^{U30)'y)/BR{B -^ K* (892)7) = 0.38 ±0.08 is found 
QQ ' in a good agreement with recent experimental data from CLEO. 

o' 
en 

^ ■ Rare radiative decays of B mesons represent an important test of the standard model of 

O ' electroweak interactions. These transitions are induced by flavour changing neutral currents 

^ ! and thus they are sensitive probes of new physics beyond the standard model. Such decays 

i-^ I are governed by one-loop (penguin) diagrams with the main contribution from a virtual top 

I I quark and a W boson. The statistics of rare radiative B decays considerably increased since 

r^' the first observation of the B -^ K*'y decay in 1993 by CLEO [|I|. This allowed a significantly 

^ . more precise determination of exclusive and inclusive branching ratios . Recently the first 
observation of rare B decay to orbitally excited tensor strange meson B -^ /Tg (1430)7 h^s 

'k>( I been reported by CLEO with a branching ratio 

BR{B -^ is:*(1430)7) = (1.66l°|^ ± 0.13) x 10-^ (1) 






as well as the ratio of exclusive branching ratios 

0.39l[;;}^ (2) 



BR{B ^ ir*(1430)7) _ , ,,+0.15 



BR{B ^ ir*(892)7) 

These new experimental data provide a challenge to the theory. Many theoretical approaches 
have been employed to predict the exclusive B -^ K* (892)7 decay rate (for a review see 
and references therein) . Less attention has been payed to rare radiative B decays to excited 
strange mesons [§-|^. Most of these theoretical approaches P,|] rely on the heavy quark 
limit both for the initial b and final s quark and the nonrelativistic quark model. However, 
the two predictions |^,H for the ratio r in Eq. (|) differ by an order of magnitude, due to 
a different treatment of the long distance effects and, as a result, a different determination 
of corresponding Isgur-Wise functions. Only the prediction of Ref. is consistent with 
data (|l|), (0). Nevertheless, it is necessary to point out that the s quark in the final K* 

1 



meson is not heavy enough, compared to the A parameter, which determines the scale 
of l/ruQ corrections in heavy quark effective theory 0. Thus the l/m^ expansion is not 
appropriate. Notwithstanding, the ideas of heavy quark expansion can be apphed to the 
exclusive B -^ K*{K2)'~f decays. From the kinematical analysis it follows that the final 
K*{K2) meson bears a large relativistic recoil momentum |A| of order of mb/2 and an 
energy of the same order. So it is possible to expand the matrix element of the effective 
Hamiltonian both in inverse powers of the b quark mass for the initial state and in inverse 
powers of the recoil momentum | A| for the final state. Such an expansion has been realized 
by us for the B — >■ i^* (892)7 decay in the framework of the relativistic quark model [Q. [| 
The obtained branching ratio for this decay was found in good agreement with experimental 
data. Here we extend this analysis to the decay B -^ Ji'g (1430)7. 

In the standard model B -^ K**'y decays are governed by the contribution of the electro- 
magnetic dipole operator O7 to the effective Hamiltonian which is obtained by integrating 
out the top quark and W boson and using the Wilson expansion ^■. 

O7 = j^sa'^'imtPR + m,PL)bF^u, Pr,l = (1 ± 75)/2. (3) 

The matrix elements of this operator between the initial B meson state and the final state 
of the orbitally excited tensor K2 meson have the following covariant decomposition 

{K;{p',e)\szKa^,b\B{p)) = ig^[e)e^^^^e*^P^k\p + p'Y . 

Mb 

{K;ip',e)\mKa,.^,b\Bip)) = g^e) [^.^iP + pX - ^If^^ip" ' P")) 

+h{e){{p' - p")k, - {p+p'),k'n,^^wz^ ^^) 

B K2 

where e^i, is a polarization tensor of the final tensor meson and k = p — p' is the four 
momentum of the emitted photon. The exclusive decay rate for the emission of a real 
photon (/c^ = 0) is determined by the single form factor g+{0) and is given by 

^(^ - "^^'^ - ^.GlmllVM.l^mn.,)KiO)§- (1 - ^)' (1 + ^) . (5) 

where Vij are the Cabibbo-Kobayashi-Maskawa matrix elements and Cj{mb) is the Wilson 
coefficient in front of the operator O7. It is convenient to consider the ratio of exclusive to 
inclusive branching ratios, for which we have 

^ BR{B -. K|(1430)7) _ 1 , Ml (l - M^jMl)' (l + Ml,jMl) 
""'^ BR{B^X,'y) r+^'Ml. {I - ml/ mlf {I + ml/ ml) ' ^' 



^It is important to note that rare radiative decays of B mesons require a completely relativistic 
treatment, because the recoil momentum of the final meson is large compared to the s quark mass. 



The recent experimental value for the inclusive decay branching ratio [|10|] 

BR{B -^ X,7) = (3.15 ± 0.35 ± 0.41) x 10"^ 

is in a good agreement with theoretical calculations. 

Now we use the relativistic quark model for the calculation of the form factor g+{0). In 
our model a meson is described by the wave function of the bound quark-antiquark state, 
which satisfies the quasipotential equation [|ll|] of the Schrodinger type [|^: 



'b\M) 



P^ 



2/i 



R, 



*Af(p) 



(2^ 



V{p,q; M)^M{q) 



where the relativistic reduced mass is 



I^R 



M' - {ml - mlf 



In the center-of-mass system the relative momentum squared on mass shell reads 



h\M) 



[M^ - {mq + mqfWM'^ - (m^ - rngf] 
4M2 



(7) 



(8) 



(9) 



The kernel Vijp, q; M) in Eq. (^ is the quasipotential operator of the quark-antiquark 
interaction. It is constructed with the help of the off-mass-shell scattering amplitude, pro- 
jected onto the positive energy states. An important role in this construction is played by 
the Lorentz-structure of the confining quark-antiquark interaction in the meson. In con- 
structing the quasipotential of the quark-antiquark interaction we have assumed that the 
effective interaction is the sum of the usual one-gluon exchange term and the mixture of 
vector and scalar linear confining potentials. The quasipotential is then defined by |1^ 



^(P, q; M) = Ug{p)uQ{-p)V{p, q; M)ug{q)uQ{-q) 
= Uqip)uQ{-p)\ -a,D^^(k)7^7^ 



-fKonf(k)r^rQ;^ + KLflk) }uqiq)uQi-q), 



(10) 



where as is the QCD coupling constant, D^i, is the gluon propagator in the Coulomb gauge 
and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 



uHp) 



e{p) + m 



/ 



1 
(jp 



X 



(11) 



V e(p) + m 
with e{p) = y/p'^ + m'^. The effective long-range vector vertex is given by 



IK 



r^(k) = 7^ + ^cr^.k'^ 



(12) 



where k is the Pauli interaction constant characterizing the nonperturbative anomalous chro- 
momagnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 



K 



^ ( 

conf \ 



'l-e){Ar + B) 



K 



conf \ 



e{Ar + B), 



(13) 



reproducing 



^conf(r) = Konf(r) + K 



V 



conf V 



Ar + B, 



(14) 



where e is the mixing coefficient. 

The quasipotential for the heavy quarkonia, expanded in v'^/c^, can be found in 
Refs. [lT3ip!4| and for heavy-hght mesons in [l^. All the parameters of our model, such as 
quark masses, parameters of the linear confining potential, mixing coefficient e and anoma- 
lous chromomagnetic quark moment k, were fixed from the analysis of heavy quarkonia 
masses |13] and radiative decays |16]. The quark masses mb = 4.88 GeV, rric = 1.55 GeV, 
rUg = 0.50 GeV, rriu^d = 0.33 GeV and the parameters of the linear potential A = 0.18 GeV^ 
and B = —0.30 GeV have the usual quark model values. In Ref. [ITTl we have considered 



the expansion of the matrix elements of weak heavy quark currents between pseudoscalar 
and vector meson ground states up to the second order in inverse powers of the heavy quark 
masses. It has been found that the general structure of the leading, ffist, and second order 
l/mq corrections in our relativistic model is in accord with the predictions of HQET. The 
heavy quark symmetry and QCD impose rigid constraints on the parameters of the long- 



range potential in our model. The analysis of the ffist order corrections |T^ fixes the value of 
the Pauli interaction constant k = —1. The same value of k was found previously from the 
fine splitting of heavy quarkonia ^Pj- states [13]. The value of the parameter mixing vector 
and scalar confining potentials e = —1 was found from the analysis of the second order 
corrections Jl^. This value is very close to the one determined from considering radiative 
decays of heavy quarkonia |T^ . 

In the quasipotential approach, the matrix element of the weak current J^ = s^k'^a^j,y{l + 
7^)6 between the states of a i? meson and an orbitally excited K2 meson has the form |jl8| 



{K*\j,mB) 



(27r)6 



^K*(p)r^(p,q)^B(q), 



(15) 



where r^(p,q) is the two-particle vertex function and '^b,k* are the meson wave functions 
projected onto the positive energy states of quarks and boosted to the moving reference 
frame. The contributions to V come from Figs. 1 and 2. The contribution F*^^) is the 
consequence of the projection onto the positive-energy states. Note that the form of the 
relativistic corrections resulting from the vertex function F^^-* explicitly depends on the 
Lorentz structure of the gg-interaction. The vertex functions look like 



F«(p, q) = w.(pi)-a^.r(l + 7')wfe(gi)(27r)35(p2 - Qs) 



and 



F!f)(p, q) = u,{p,)u,{p,)\Ua,,,K{l + 7i^) ,t\ 1_^'] , 7i°V(p2 - q2 



+V(P2 - q2: 



es{k[) + es{qi) 



h) + efe 
liicrifj,^ku{l + 7i) \ub{qi)uq{q2), 



(16) 



(17) 



where ki = pi - A; k'^ = qi + A; A = pxj - Pb, 



A^-\p) 



e{p) - (m7° + 7°(7p)) 



It is important to note that the wave functions entering the weak current matrix ele- 
ment (^) cannot be both in the rest frame. In the B meson rest frame, the K2 meson is 
moving with the recoil momentum A. The wave function of the moving K2 meson "^k* /^ is 
connected with the K2 wave function in the rest frame "^k* = "^k* by the transformation 



*i.| a(p) = Dy\RfjDl/'{RYj<i^K.^o{p) 



(18) 



where R^ is the Wigner rotation, La is the Lorentz boost from the meson rest frame to a 
moving one. The wave functions of B and K2 mesons at rest were calculated by numerical 
solution of the quasipotential equation (|^. 

We substitute the vertex functions T^^^ and F*^^) given by Eqs. (^) and {\n\) in the decay 
matrix element (|T3p and take into account the wave function transformation (^H|). The 
resulting structure of this matrix element is rather complicated, because it is necessary to 
integrate both over d^p and d^q. The 6 function in expression (|T6|) permits us to perform one 
of these integrations and thus this contribution can be easily calculated. The calculation of 
the vertex function F^^^ contribution is more difficult. Here, instead of a 5 function, we have 
a complicated structure, containing the qq interaction potential in the meson. However, we 
can expand this contribution in the inverse powers of the heavy b quark mass and large 
recoil momentum | A| ~ mii/2 of the final K** meson. Such an expansion is carried out up 
to the second order. ^ Then we use the quasipotential equation in order to perform one of 
the integrations in the current matrix element. As a result we get for the form factor g+{0) 
the following expression with k, = —1 



940) = 9^^\0)U1 - e)^f "'(0) + 5(7f ^(0) 



(19) 
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1 + 



M, 
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Mb - Ek' 



+ 



p 



p 



eq{p)+mq es{p + A)+ms 
P' 



es{p + A) +ms [es{p + A) + ms][ebip) + mh\_ 



>^Bip), 



(20) 



^This means that in expressions for g_^_' (0) and g]^_ (0) we neglect terms proportional to the 
third order product of small binding energies and ratios p^/eg(A), p^/ej^(A) as well as higher order 
terms. 
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pA2 2[e,(A) + m,] eq{p) + m, 

(Ms + Mk*)[Mb - efe(p) - eg(p)] + {Ek* + M^j) 
2e. 



E^. + M^. 



'B -C/KJ + Mk* 






- e 


i(p + 


Eu'* 


''"• All 


(p 


+ F, 


2eg 


v,...^i 




(21) 
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eg(A) + mg 
2e,(A) 
eg(A) -ms 
2e,(A)[e,(A)+m, 
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K* 
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(p + 
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Mb - ^i^* 
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<(P + 
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.(p + 



Ek-^ + M^* 



Mb + Mk' - eb{p) - eg{p) 



>i^B{p) 



(22) 



where the superscripts "(1)" and "(2)" correspond to contributions coming from Figs. 1 and 
2, S and V mean the scalar and vector potentials in Eq. ([13D , 4'k*,b are radial parts of the 
wave functions. The recoil momentum and the energy of the K2 meson are given by 



M|-Mi,_ 



2M, 



E 



M% + Ml 



K* 



K* 



B 



2Mf 



(23) 



We can check the consistency of our resulting formulas by taking the formal limit of b 
and s quark masses going to infi nity. P| In this limit according to the heavy quark effective 
theory [|l^ the function ^p = 2 J Mb Mk* g+/{ Mb + Mk^) should coincide with the Isgur- 
Wise function r for the semileptonic B decay to the orbitally excited tensor D meson, 
B -^ Dj^i/. Such semileptonic decays have been considered by us in Ref. [0. It is easy to 
verify that the equality of ^p and r is satisfied in our model if we also use the expansion 
in {w — l)/{w + 1) {w is a scalar product of four-velocities of the initial and final mesons), 
which is small for the B — > D^^ev decay [^]. Using Eq. ([T9|) to calculate the ratio of the 
form factor 5'+(0) in the infinitely heavy h and s quark limit to the same form factor in the 
leading order of expansions in inverse powers of the heavy h quark mass and large recoil 
momentum | A| we find that it is equal to Mb/ JM^ + Mj^* ^ 0.965. The corresponding 

ratio of form factors of the exclusive rare radiative B decay to the vector K* meson -Fi(O) (see 
Eq. (23) of Ref. [|) is equal to Mb/ J Ml + M|. ^ 0.986. Therefore we conclude that the 



^As it was noted above such limit is justified only for the h quark. 



ratios of form factors 5'+(0)/-Fi(0) in the leading order of these expansions differ by factor 



M]^ + M]^t / J M^ + Mj^* ^ 0.98. This is the consequence of the relativistic dynamics 
leading to the effective expansion in inverse powers of the s quark energy es{p + A) = 



p + A)2 + ml, which is large in one case due to the large s quark mass and in the other 
one due to the large recoil momentum A. As a result both expansions give similar final 
expressions in the leading order. Thus we can expect that the ratio r from (^ in our 
calculations should be close to the one found in the infinitely heavy s quark limit 0. 

The results of numerical calculations using formulas (||), (P), (p!9D -(p2|) for e = —1 are 
given in Table |. There we also show our previous predictions for the B —>■ K*'~f decay [Q . 
Our results are confronted with other theoretical calculations [§-0] and recent experimental 
data 0. We find a good agreement of our predictions for decay rates with the experiment 
and estimates of Ref. [Q. Other theoretical calculations substantially disagree with data 
either for B -^ K*'~f Q or for B -^ -R'It B decay rates. As a result our predictions and 
those of Ref. |^ for the ratio r from (Q) are well consistent with experiment, while the r 
estimates of [§] and |^ are several times larger than the experimental value (see Table |). As 
it was argued above, it is not accidental that r values in our and Ref. approaches are close. 
The agreement of both predictions for branching ratios could be explained by some specific 
cancellation of finite s quark mass effects and relativistic corrections which were neglected 
in Ref. 0. Though our numerical results agree with Ref. 0], we believe that our analysis is 
more consistent and reliable. We do not use the ill-defined limit rris —>■ oo, and our quark 
model consistently takes into account some important relativistic effects, for example, the 
Lorentz transformation of the wave function of the final K** meson (see Eq. (|18D). Such a 
transformation turns out to be very important, especially for B decays to orbitally excited 
mesons [^. The large value of the recoil momentum | A| ~ mb/2 makes relativistic effects 
to play a significant role. On the other hand this fact simplifies our analysis since it allows 
to make an expansion both in inverse powers of the large b quark mass and in the large 
recoil momentum. 
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TABLES 



TABLE I 
for branching ratios and their ratios Rk 
values for the S 



Our results in comparison with other theoretical predictions and experimental data 

BRjE^K'-y) jj _ BR{B^K*'y) ^_ BR{B^K*'y) 



BR{B-*Xs'y) 

K*'y decay are taken from Ref. [8]. 



P ^ BR{B- 
^K* — BR{B~ 



-Xs-i) 



BR{B-^K*y) ■ 



Our 



Value 



our 



Ref. g Ref. [| 



Ref. § 



Exp. § 



4:.55tof8 ± 0-34" 



BR{B^K*j)xW^ 4.5 ±L5 



Rk* (%) 
BR{B -^ KIj) X 10^ 

Rk; (%) 



r 



15 ±3 

L7±0.6 

5.7 ±L2 

0.38 ±0.08 



L35 

4.5 
L8 
6.0 
1.3 



1.4-4.9 

3.5-12.2 

6.9-14.8 

17.3-37.1 

3.0-4.9 



4.71 ±1.79 

16.8 ±6.4 
1.73 ±0.80 

6.2 ±2.9 
0.37 ±0.10 



3.76' 



-0.89 
-0.83 



± 0.28" 



1.66to:53±0.13 



0.39 



-0.15 
-0.13 



" B^ ^ K*^j 
^ B+ ^ K* 



7 



7 



FIGURES 



^ ^ 


\ -"^ 


B 


K*2 


q 


q 



FIG. 1. Lowest order vertex function F^^' corresponding to Eq. (ITI 
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FIG. 2. Vertex function F"^ corresponding to Eq. (p!7[). Dashed lines represent the 
effective potential V in Eq. (10). Bold lines denote the negative-energy part of the 
quark propagator. 



10 



